Bends in Optical Dielectric Guides 

By E. A. J. MARCATILI 

(Manuscript received March 3, 1969) 

Light transmission through a curved dielectric rod of rectangular cross 
section embedded in different dielectrics is analyzed in closed, though 
approximate form. We distinguish three ranges: 

(i) Small cross section guides such as a thin glass ribbon surrounded 
by air — Making its width 1 percent of the wavelength, most of the power 
travels outside of the glass; the attenuation coefficient of the guide is two 
orders of magnitude smaller than that of glass, and the radius of curvature 
that doubles the straight guide loss is around 10,000\. 

(ii) Medium cross section guide for integration optics — It is only a few 
microns on the side and capable of guiding a single mode either in low loss 
bends with short radii of curvature or in a high Q closed loop useful for filters. 
Q's of the order of 10 s are theoretically achievable in loops with radii ranging 
from 0.04 1° 1 mm, if the percentage refractive index difference between guide 
and surrounding dielectric lies between 0.1 and 0.01. 

{Hi) Large cross section guides — They are multimode and are used in 
fiber optics. Conversion to higher order modes are found more significant than 
radiation loss resulting from curvature. 

I. INTRODUCTION 

A dielectric rod, embedded in one or more dielectrics of lower re- 
fractive index, is the basic ingredient of three types of optical wave- 
guide which differ only in their relative dimensions and consequently 
in their guiding properties. 

The first is a small cross section guide which supports only the fund- 
amental mode; most of the power travels in a lower loss external 
medium. Thus, the attenuation of the mode is smaller than if all the 
power flowed through the higher loss internal medium. Tiny rods, thin 
ribbons, or films made of glass or other substances embedded in either 
air or low loss liquids are typical examples. 1-3 

The second is a medium size guide capable of supporting only a few 
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modes; most of the power travels in the internal medium. Such a guide, 
(Fig. 1 of Ref. 10) has been proposed as the building block of passive 
and active components for integrated optical circuitry. 4 " 6 Lasers, 
modulators, directional couplers, and filters are some of the many 
devices which could be built in a single substrate utilizing the high 
precision techniques available from integrated circuitry; consequently 
they would be compact, mechanically stable, and reproducible. 

The third, a large size guide (clad fiber) which can support many 
modes, is used typically in fiber optics. 7 

These basic guides, having round or rectangular cross section and 
straight axis, have been studied both analytically and through com- 
puter calculations. 8 " 13 Also the directional coupler (Fig. 2 of Ref. 10) 
obtained by running two guides of rectangular or circular cross sections 
parallel to each other, has been analyzed. 10 - 12 - 14 

To my knowledge, though, little is known quantitatively about the 
ability of any of the three types of guides to negotiate bends, or about 
the radiation losses in loops, such as the one depicted in Fig. 1 as part 
of a channel dropping filter. This paper should supply such informa- 
tion. 

In Section II the boundary value problem is discussed, and the 
fundamental modes of each polarization are described. Section III con- 
tains a discussion of the results and numerical examples. Conclusions 
are drawn in Section IV and all the mathematical derivations are 
exiled to the appendix. 

II. FORMULATION OF THE BOUNDARY VALUE PROBLEM 

Figure 2 depicts, in perspective, the basic geometry of the curved 
guide with radius of curvature R. The cross section is a rectangle whose 
sides are a and b. The refractive index of the guide is n x , and the re- 
fractive indices around the guide are no, n 3 , n 4 , and n 5 , all of which are 
smaller than rii . Furthermore, for reasons which become apparent later, 
we do not specify the refractive indices in the four shaded areas. 

This boundary value problem is solved in closed, though approxi- 
mate form in the appendix, by introducing the same simplification 
used in solving the problem of transmission in the straight guide. 10 
That simplification arises from solving Maxwell's equations only for 
guide dimensions such that a small percentage of the total power flows 
through the shaded areas and consequently a negligible error is ex- 
pected if one does not match properly the fields along their edges. 

Two types of hybrid modes propagate through this curved guide; 
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Fig. 1 — Channel dropping filter (ring type). 

each one has six field components. But since some of the refractive 
indices n 2 , n 3 , n 4 , and n 5 are chosen close to n r , guidance occurs 
through total internal reflection only when the plane wavelets that make 
a mode impinge on the interfaces at grazing angles. Consequently, the 
only large field components are perpendicular to the curved z axis 
(Fig. 2). The modes are then of the tem kind and we group them in two 
families, E' va and E V VQ . The main field components of the members of the 
first family are E x and H v , while those of the second are E v and H x . 




Fig. 2 — Curved dielectric guide. 
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Virtually every one of these components varies sinusoidally along 
x and y within the guiding medium 1 and decays exponentially in the 
surrounding media 2, 3, 4, and 5 (Fig. 2). The subindices p and q 
represent the number of extrema of each field component in the x and 
y directions, respectively. The field configurations of some members 
of the two families in straight guides are depicted in Fig. 5 of Ref. 10; 
section 2.1 describes the influence of a finite radius of curvature on 
those field configurations. 

General expressions for the different phase and propagation constants 
in each medium of the curved guide are calculated in the appendix, for 
arbitrary modes and for n 2 ^ n 3 -^ n 4 5^ n s . In the text, we consider only 
the fundamental modes of each family E x n and E" n ; furthermore, we 
choose 

n 3 = n 5 (1) 

and leave ru, and n A arbitrary. This choice of refractive indexes en- 
compasses the most interesting cases. 

2.1 E x n Mode 

We first study the E X 1X mode. As we said before, the main compo- 
nents are E x along the x direction and H v along y. Both components 
have a single maximum located within medium 1 and drop sinusoidally 
toward the edge of it. Outside of the medium, the decay is exponential. 

The axial propagation constant is according to equation (47) 

K = {k\ -kl- #)», (2) 

where fci = /c% = (27r/A)% and A is the free space wavelength, k w is 
the propagation constant along x in media 1, 2, and 4, and k v is the 
propagation constant along y in media 1, 3, and 5. This means that 
the electrical width of media 1, 2, and 4 is the same and equal to k x a, 
and the electrical height of 1, 3, and 5 is also the same and equal k v b. 
The transverse propagation constant k y is independent of the radius 
of curvature R and can be found from the transcendental equation 
(37) 



k v b = t — tan" 
in which 



\k u A 2 ) 1 \ 



— tan" 



_ V/CyA^ 



(3) 



4 2[n\ - nty 
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If the height of the guide b is selected so large that 

Ai \ Ai « 1, (5) 

only a small percentage of the power carried by the mode travels in 
media 2 and 4; and equation (3) can be solved approximately, 
yielding 

According to equation (49), the other transverse propagation con- 
stant 

fc --4 1 + ^- i tf (6) 

is valid if 

c 



a c R « 1. 



(7) 



The first term in equation (6), k z0 , is the propagation constant in the 
x direction of the guide without curvature; the second and third terms, 
which according to equation (7) must be small, are perturbations related 
to the change of field profile and to radiation loss, both of which are 
introduced by the curvature. More precisely, a e is the attenuation 
coefficient of the curved guide, a c R is the attenuation per radian, 
that is the attenuation in a length of guide equal to R, and c is a con- 
version loss coefficient such that, at a junction between a straight and a 
curved section of the same guide, c measures the power that the funda- 
mental mode in the straight section would couple to modes higher than 
the fundamental in the curved section. The fact that equation (6) is 
valid if c <$C 1 requires the radius of curvature R to be so large that the 
field profiles of the fundamental modes in the straight and curved guides 
are quite similar. Later in this section we consider formulas applicable 
when c £g 1. 

The axial propagation constant, k l0 , of the straight guide is related 
to k x0 and k v by the expression 

*.«, = (ti ~ fc*o - &; (8) 

and k x0 is the solution of the transcendental equation (55) 
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k in a = 



^-IteJ- 1 ! 



The length 



(9) 



A = 



2(n? - n\) 



»\l 



(10) 



is used as a normalizing dimension. What does it measure? If one 
assumes b = », the guide becomes a slab of width a. If a S A, only 
the fundamental mode is guided; if a > A, the slab is multimode. 

Figure 3 is a graph of the electrical width, k x0 a, of the straight guide 
as a function of a/ A. The solid curve is the solution of equation (9) 
assuming n x /n 3 = 1.5, while the dotted one is the solution for n,/n 3 =1. 
For thin guides, a/ A « 1, the electrical width is proportional to a; for 
thick guides, a/ A » 1, the electrical width goes asymptotically to w. 

The attenuation per radian a c R and the conversion coefficient c, 
obtained from equations (50) and (51) with n 3 = n s are 



M = 



1 - 



f#-* 



flslOxod 



1 - 



k T(i A 



IT 



(R exp S — jr 



2c 



1 -'¥a i+ «*.j 



-(-:i)(¥r-f [-(¥)'] 



-J 



(11) 



and 



c = 



1 Iwa 

2fc l0 a \i4 



1_ 

(R 



(12) 



where 



2^ ALA^Yfl 

^ " k 2 t0 A 3 ~ 2 fc!o V 1 n 2 J K - 



(13) 



The solid curves in Figs. 4 and 5 are graphs of the function 



2 \* 



a c R[ 1 - 



(which is proportional to the attenuation per radian) as a function of a/ A 
using (Rasa parameter. In Fig. 4, we further assume that 



^= 1 + A 
n 3 
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NORMALIZED WIDTH -j- = ^ 



Fig. 3 — Guide's electrical width. Solid line is for E^ mode with njnt = 1.5; 
dashed line is for S", mode with ni/n 3 arbitrary, and for E x n mode with n\ ^ n 5 . 



and 



in Fig. 5, 



A« 1; 



» = 1.5. 

n 3 

In the same figures each dashed line is a curve of constant conversion 
loss c. Since the calculations are valid for c « 1, we believe the solid 
curves are reliable to the left of the dotted curve c = 0.3 and grow 
progressively in error to the right of it. 

To extend the use of this graph to arbitrarily large values of a/ A, we 
calculate the loss per radian, equation (63), when a/ A » 1 andc^ 1. It is 



a r R = 



nil 



l-l* 



-i 



exp 



-![ 



^-<l 






(14) 



the dotted lines in Figs. 4 and 5 represent this loss. The reader can 
smoothly extend the solid curves to the right of the dashed line, c = 
0.3, so that they become asymptotic to the dotted lines. Thus, the 
whole range of guide width a from to R has been covered. 
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Fig. 4 — Attenuation per radian for £?n s and En" modes if ni/rii = 1 -f- A 
and A « 1. 

To understand why these curves of constant R become asymptotic 
for a/A 3> 1, we have drawn in Fig. 6a a curved guide with a certain 
R; its width a is very large compared with A. Also the amplitudes of 
the field components E x and H y are plotted as functions of x and y. 

Along x the field inside the guide behaves virtually as the Bessel 
function J v [k x (R + x)] where v is a very large number and outside 
of the guide decays exponentially. This guide has some radiation loss 
per radian. 

Now, suppose that we start shrinking a without changing R. Since 
the field at x = —a is very small, the radiation loss remains constant 
until a is made so short that the field at x = and x = — a are com- 
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parable (Fig. 6b) . The field inside the guide varies almost sinusoidally, 
while outside decays exponentially and the attenuation per radian 
increases. If a is reduced even further (Fig. 6c) most of the power 
travels outside of the guide, and the loss increases even more. The 
field configuration along y is practically the same in the three cases 
(Fig. 6). 

For resonant loops, such as the filter in Fig. 1, the intrinsic Q re- 
sulting from curvature radiation is more interesting than the attenua- 
tion a c . They are related by the expression 
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Fig. 5 — Attenuation per radian for S n mode if ni/n 3 = 1.5. 
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Fig. 6 — Field distribution as a function of guide width a with (a) a/A » 1, 
(b) a/A ^ 1, and (c) a/A « 1. 

This function is plotted in Fig. 7, assuming 



l -i = 1 + A 



71 a 



and 

and in Fig 8, assuming 



A«l 



^ = 1.5, 
n 3 



BEND.S IN GUIDES 



2113 



using as before the normalized guide width a/ A as variable and (R as 
parameter. As in Figs. 4 and 5, the reader can easily match the solid and 
dotted curves. Further discussion of these curves is reserved for Section 
III. 

The field components in media 2, 3, 4, and 5 decay almost exponen- 
tially away from the guiding rod, and the distances ??._,, 774, £3, and £ 5 
over which the fields decrease by 1/e are 
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Fig. 7 — Intrinsic Q for En" and En" modes if ni/n 3 = 1 + A and A « 1. 
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2T^' 



(16) 
(17) 



2.2 £?, Mode 

We now consider the 2SJ, mode. The main components are £"„ and # x ; 
they are qualitatively quite similar to components of the E x n mode, 
rotated 90°. 

The propagation constant k t is still given by equation (2); but now k y 




Fig. 8 — Intrinsic Q for En* mode if ni/rij = 1.5. 
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is the solution of 



k„b = 



tan -5 



n\ \_\k y 



- 1 



— tan 



2 



_ \'l-*y- / l4 



- 1 



-i 



(18) 

The equivalent formula of any of those between equation (7) and (17) 
can be derived from that formula by substituting the ratio of refractive 
indexes by unity, but leaving them unchanged wherever they are 
subtracted from unity. For example, equation (11) becomes 



«■*-£ 



!_ Si 
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(KoO)'' 
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TT(l 
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ICj-qA 



(R exp S — 



(J{ 



1 - 



k ia A 

■K 



1 + 



2c_ 
ak x J J 



1-21- 



n 3 



k x0 A\ 



>ib-ffl. 



(19) 



while c and (R given by equations (12) and (13) remain unchanged. 

Figure 9 is a graph of the function a c R[l — (n 8 /ni) 2 ]*, valid for any 
ratio rii/na ■ In particular, for njn 3 = 1 + A and A « 1, equations (19) 
and (11) become the same, and consequently these curves coincide with 
those in Fig. 4. This means that for w, ^ n 3 , the E x n and E^ modes have 
the same loss. 

Figure 10 is a graph of the intrinsic Q of a loop operating in the E v tl 
mode which can be derived from equations (15) and (19). As before, in a 
resonant loop with n x ln z = 1 + A and A « 1, the E\ x or E" n modes have 
the same Q's. 



III. DISCUSSION AND EXAMPLES 

The attenuation per radian of any dielectric guide of rectangular 
cross section and the Q c resulting from curvature are strongly de- 
pendent on the radius of curvature. With the help of equation (17), 
the attenuation per radian equation (11) can be written 



a c R = MR exp I - 



1 \IR 



(20) 



6V | £ a . 

where M is independent of R, \. is the guided wavelength along z, and 
£j is the length over which the field in medium 3 decays by 1/e. Ac- 
cording to Fig. 11, the function 

1 \ 2 ,R 



Rexp[-— 2 



kz 
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becomes negligibly small, and consequently the attenuation per radian 
becomes negligibly small when 



R > 



24tt 2 1 £ 3 



(21) 



This simple criterion is developed further in Ref. 15. 

We are interested, though, in a more detailed description of trans- 
mission through a bent dielectric guide. Given a guide with a certain 
radius of curvature (that is, given R and a/ A), in general the loss per 
radianof the E\ x mode is much larger than that of the E\ x mode (compare, 
for example, Figs. 5 and 9 for n x jn z = 1.5). That difference becomes 
negligible if njn 3 — 1 <$C 1. 
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Fig. 9 — Attenuation per radian for Eu v mode and ni/n 3 > 1. 
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Fig. 10 — Intrinsic Q for En* mode and ni/n 3 > 1. 

Let us consider separately the three types of guide: thin, medium and 
large. 

3. l Thin or Low Loss Guides* 

In thin guides the width a is so small that 



* Low loss for straight guide. 
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a 2a(n] - n§* 



A 



A 



«1. 



(22) 



The height 6 of the guide must be large so that only a little part of 
the power travels in the shaded areas of Fig. 2. Assuming that the 
guiding rod dielectric is lossy, its refractive index is 



h = n ( x + £[) ' 



(23) 



where n is real and a is the attenuation constant of a plane wave in 
that medium. 

Substituting equations (22) and (23) in equations (2), (11), and 
(12), we obtain 

k z = k in + ia, + i<*c ■ (24) 

The first term 



k, = (kt - k 2 u )K 



l + l 



fc,al 1 — 



1 + 



1 



n 






for E' n mode 



for £"11 mode 



(25) 



is the phase constant. Since most of the power travels in the external 
medium, its value for either mode is close to /cn 3 . The conversion loss 
term c is negligible. 

The imaginary part of equation (24) is the attenuation constant, 
and is made of two terms. The first term 



a 72 2 in , w 



— I for E x tl mode 
1 for E" u mode 



(26) 



£-efc) 




R HR 4-Rq 

Fig. 11 — Plot of R/Rq exp ( —R/R ) and tangent at inflection point. 
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is the attenuation that each mode would have if the guide were 
straight. 10 The second term 



kW 


(iH* 




p^m-vb-m) 




(for E x n mode) 




-{-WS-iJfr-iOT 




(for Elt mode) 



(27) 



is the attenuation resulting from the radiation introduced by the curva- 
ture. The E" n mode is more tightly bound to the guiding rod and conse- 
quently has more straight loss and less curvature loss than the E x u mode. 
From equations (26) and (27), the radius of curvature R d that doubles 
the straight guide loss is 



k 3 



an 



2a, n : . 



n ->-\ 



'!>-§( 



?j iog 



. 



\kn l\ 

[knl (n 
log \&k \n 



2.-1 



- 1 



(for E x n mode) . 



(for E" u mode). 



(28) 



Example 1 : Consider a thin ribbon guide made of glass surrounded by 
air and assume that n = 1.5, n 3 = 1, a = 0.1 nepers per m, and b = °o . 
From equations (26) and (28) we calculate the values in Table I. 

It is doubly advantageous to use the E\ x mode rather than the E\ x 
because (i) the thickness required for equal radiation loss and straight 
guide loss is roughly (n/n 3 ) 3 times larger, and (n) R d is about (n/n 3 )* 
times smaller. 

If the height b of the ribbon is finite, k v /kn 3 is no longer zero and the 
radii are, according to equation (28), [1 — i(fc»/fcO*T 1 times longer than 
those in Table I. 

3.2 Medium Size Guide for Integrated Optical Circuitry 

It is likely that guides for integrated optical circuitry will be possible to 
fabricate only with n, ^ n 3 . The radiation loss per radian and the Q e of 
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Table I — Values Calculated from Equations (26) and (28) 





Euv Mode 


Eu* 


Mode 


(nepers /m) 


a Rd 
X X 


a 
X 


Rd 
X 


0.01 

0.001 

0.0001 


0.05 1.9 X 10 3 
0.016 6.2 X 10 4 
0.005 2 X 10 8 


0.17 

0.055 

0.017 


6.3 X 10' 
2 X 10* 
6.5 X 10 6 



loops made with these guides can be obtained from Figs. 4 and 7, con- 
sidering abscissas around a/ A = 1. For both modes, 2?,, and E' u , most 
of the power travels within the guiding rod.* 

In general, the losses are very sensitive to the radius of curvature. 
They are also sensitive to the guide's width to the left of the dashed 
curve c = 0.5, but fairly insensitive to the right of it. 

Example 2: Let us design a guide : 

(i) The attenuation per radian resulting from radiation loss is 

a c R = 0.01 nepers = 0.087 dB. 

(ii) Its width a is the maximum compatible with single mode guid- 
ance in the infinitely high slab, that is 

d za, ,2 2\* t 

(Hi) We assume b = <x> and n 3 = TOi(l — A), where A « 1 andn, = 
1.5. 
From Fig. 4 we derive the guide dimensions for different values 

of A: 



0.1 

0.01 

0.001 



a 

x 

0.745 

2.36 

7.45 



R 
X 

30 

1,060 
37,000 



Unless A is 0.01 or larger, the radius of curvature R becomes un- 
comfortably large for integrated optical circuitry. Furthermore, if b 
is finite, k y is no longer zero, and the radii become [1 — (/^/fc.-,) 2 ] -1 
times larger than those in the table above. 



* This is not true if b/B 2 « 1. Then k,o must be calculated from equation (8). 
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Example 3: We design a resonant loop (Fig. 1) such that its Q c re- 
sulting from radiation is equal to the Q resulting from transmission 
loss in typical glass (w, = 1.5, a = 0.1 neper/m at A. = lfi) ; that is, 

Q = Q c = 5 X 10 7 . 



Furthermore, let us assume as in Example 2 that a/A = 
n,(l — A), and b = <x> . With the help of Fig. 7 we derive 



1, n a = 



A 

0.1 

0.01 

0.001 



a 
X 

0.745 

2.36 

7.45 



R 
X 

57 

1 , 550 

42,000 



Again, unless A is larger than 0.01, the radius of curvature becomes 
unwieldily large for integrated optical circuitry. 

Instead of using a loop as the resonant circuit of Fig. 1, it is pos- 
sible to make a = R, and the loop becomes a pillbox (Fig. 12). This 
structure may be simpler to fabricate. For this case, also from Fig. 4, 
using the refractive indices of the previous example, we obtain 



A 

0.1 

0.01 
0.001 



R 
X 

42 

1,170 

32,000 



The pillbox resonator requires a 30 percent shorter radius than the ring 
resonator. As before, if b is finite, the radii are [1 — (k y /k 3 ) 2 ]~ 1 times 
longer than those in the last two tables. 



PILLBOX RESONANT 
AT f, 



/v 3 




Fig. 12 — Channel dropping filter (pillbox type). 
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3.3 Large Guides for Fiber Optics 

The large guide is multimode, a/A » 1, and the radius for small 
mode conversion is derived from equations (11) and (12), making 
k x0 a = it and k e0 = 2ir?ii/A. Then 

.3 

2 a 

For a power conversion c 2 = 0.01, and Wi = 1.5, we have 

a B 

x x 

5 8,900 

10 71,000 

The conversion loss is many orders of magnitude larger than the loss 
radiated by the fundamental mode because of the curvature. Radia- 
tion loss of higher order modes can be found in equations (51) and 
(63). 

In general, clad fibers are of circular cross section; consequently 
our calculations do not strictly apply. Nevertheless, a guide of circular 
cross section and another of equal area but square cross section must 
have quite comparable attenuation per radian unless mode degen- 
eracy occurs, but this is quite unlikely. 

Though we have been talking throughout of light guides, it is 
obvious that all the calculations are equally applicable to microwave 
guides. 

IV. CONCLUSIONS 

Relations between radiation losses resulting from curvature, geom- 
etry, and electric characteristics of the bent dielectric guide are sum- 
marized in Figs. 4, 5, and 7 through 10 and they are discussed and 
exemplified in Section III. 

The main qualitative results are that for a given radius of curva- 
ture R, the radiation loss can be reduced 

(i) by increasing the difference between the refractive index ni of 
the guide and those of the media toward the outside, n 3 , and inside, 
715, of the curved guide axis (Fig. 2) ; 

(w) by increasing the guide width a. Nevertheless, once a is bigger 
than 
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(where A is the free space wavelength), there is little reduction of the 
loss ; 

(Hi) by choosing the height of the guide large enough to confine 
the fields as much as possible within the guide in the direction normal 
to the plane of curvature. 

In general, the radiation losses are small if 

R > ^lk! 
« > x , 

where £ 3 is the length over which the field decays by 1/e in medium 3 
(Fig. 2). 

Thin ribbons of glass, surrounded by air and oriented as in Fig. 6c, 
operate better with the electric field perpendicular to the ribbon's 
plane. Choosing the thickness a = 0.055A, the attenuation of the 
straight guide is 1 percent of the attenuation in glass, and the radius 
of curvature which doubles that low attenuation is 20,000A. 

The dielectric guide for integrated optical circuitry seems suitable to 
negotiate bends and to make resonant loops of small radii of curvature 
and small radiation losses. For example, for 

n x = 1.5 
a = zr-. (single mode guide) 



2»,| 1 - -i 



a 1 percent attenuation (0.087 dB) resulting from radiation in a length 
of guide equal to R is achieved with the following values 

. _ nj a R 

71-1 X \ 

0.1 0.745 30 

0.01 2.36 1060 

0.001 7.45 37000 

The smaller ?ii — ?i 3 , the larger the radius of curvature. For A = 0.63/x, 
if one wants to keep R below 1 mm, the difference between the internal 
and external refractive indices must be larger than 0.01. 

Large cross section dielectric guides capable of supporting many 
modes are far more sensitive to mode conversions than to radiation 
losses. For the fundamental mode, the power conversion loss at the 
junction between a straight and a curved section of a multimode 
guide is 
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a \2 

, 2 a I 

For ri! = 1.5, a = 6.3/x, and A = 0.63/*, the radius of curvature R that 
produces a power conversion c 2 of 0.01 is 45 mm. The radiation loss in 
a length of guide equal to R is many orders of magnitude below 0.01. 

APPENDIX 

Field Analysis of the Curved Guide 

Figure 2 shows the geometry and dielectric distribution of the curved 
guide. In this appendix two families of modes are found, E x va and E v vq ; 
each is studied separately. 

a.i E x pq Modes: Polarization Along x 

The field components in each region should be written as integral 
expressions, but, as discussed in Section II, the power propagating 
through the shaded areas is neglected, and the field matching is per- 
formed only along the sides of region 1. Consequently, those field 
components do not need to be so general. As a matter of fact, the 
simplest field components in the rath of the five areas are 16 

i 



H x „ 



1 d'H u 



K, - k; m dx dy 



H um = e 



-i»«+i 



Hzm — 



\M l J.[(k* l - k 2 vl ) h (R + *) + to] cos (k vl y + 0,) for m = i 

M 2 J, [\k\ - k 2 v2 y(R + x) + to ] exp f =Fik ya y'\ for m = * 

J 

M 3 Hl 2) [(k 2 3 - k 2 yZ )KR + x)] cos (k u3 y + fi 3 ) for m = 3 

[M t J,[(kl - kl B )*(R + »)] cos (k vl y + Q fi ) for m = 5 

i v dH um 

kl- k 2 yn R + x dy ' 



**■ ~ k 2 m -kl m (R + x) a " m ' 

E vm = 0, 

*- fci-fcL dx ' ^ } 
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in which M m is the amplitude of the field in the mth medium; ^ m and £2 m 
are constants that locate the field maxima in region m; « is the angular 
frequency; en 2 and /x, the permittivity and permeability of each medium, 
are related by k 2 = k 2 n 2 = ofenn 2 ; k um is the propagation constant 
along y in medium m; and J, and H ( , 2 * are Bessel and Hankel functions, 
respectively. 

Strictly speaking, the H y component in media 1, 2, and 4 should 
be written as a sum of Bessel functions of the first and second kind, 
but later on they are approximated by circular functions; therefore, 
we do not make any mistake using only the Bessel function of the 
first kind with an arbitrary phase constant in the argument. 

We consider only guide geometries for which the guide wavelengths 
measured in the x and y directions in medium 1 are large compared 
with the wavelength measured in the z direction. This means that (i) 

dJ t«k' ™ 

and, as a consequence, the field component H x \ is very small com- 
pared with H z and is neglected; (u) the propagating modes are bas- 
ically of the tem type. 

In order to match the remaining components along the boundaries 
of medium 1, the field components in media 1, 2, and 4 must have the 
same dependence along x, while the field components in media 1, 3, 
and 5 must have the same dependence along y. Therefore 

K ul = AC„3 = /C„ 5 = K v , (31) 

^i = & = & = ^, and a x = fi 3 = fis = n. (33) 

Furthermore, the field matching yields the following four equations 
from which two characteristic equations will be derived 

b . „\ . A*„a I , b „\ . k ui 



tan l*» 5 + Oj - * IT i tan \ k » 2~ n ) = { k ^ 34) 

and ^ = ^^1 (35) 

JAPia) Pis J ,{pn) 





J,{p„) 

j:( p > 3 ) 


Pa H l , 2) (p 3 ) 

p, 3 Hi "(p 3 ) 


where 








Pta 


= R(k'i - 


kl)> 


+ *, 


Pi 


= R(kl - 


kl)K 


and 



Pl5 = (R - a)(k 2 - k 2 v f + 4, 
P5 = (R - a)(k 2 5 - kl)*- 



(36) 
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Similar to what happens with the straight guide, equations (34) 
and (35) are the boundary conditions of two independent problems 
far simpler than the one depicted in Fig. 2. Thus, for a dielectric slab 
infinite in the x and z directions and with dimensions and refractive 
indices as depicted in Fig. 13a, the boundary conditions for modes 
with no E y component coincide with equation (34). Similarly, for a 
bent slab infinite in the y direction as shown in Fig. 13b, the boundary 
conditions for modes with a negligible H T component coincide with 
equation (35). 

The elimination of fi between the two expressions of equation (34) 
yields the characteristic equation for the plane slab 10 



k u b = qir — tan" 



in which 



1 



A?k< 



- 1 



— tan ' f 



1 



(aiJ - ']' ' 



(37) 



A, = 



2\n 2 x — nly 



(38) 



the tan -1 functions are to be taken in the first quadrant, and the 
arbitrary integer q is the order of the mode, that is, the number of 
extrema of each field component within the guiding rod in the y 
direction. 

The transcendental equation (37) has an approximate closed form 
solution already found in Reff. 10 



*. » ¥■ 1 



cpr 

b 



A 2 + A, 



+ 



(39) 




(a) (b) 

Fig. 13 — Guiding dielectric slabs. 
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(40) 



and consequently the parenthesis is close to unity. 

The field components in media 2 and 4 decay exponentially by 1/e 
in lengths r/ 2 and 774, which are deduced from equation (32) to be 



v-2 = 



IC U 2 I "»1 "-2 "'j; ) 

4 I \ 4 / 



(41) 



Let us consider the solution of the characteristic equation of the bent 
slab (Fig. 13b) . For guided modes, both the arguments and the order 
of the Bessel and Hankel functions involved in equation (35) are large 
compared with unity, and consequently they can be replaced by their 
Watson's first term approximations, 17 



J.(p) = 



Up* - , 2 ) 5 J 



exp 



/ 2 2\i 

{v — py 



(p 2 - S) 



3k 

2\l 



+ 



Y.(p) = - 



_tt(p — v y _ 
These expressions are valid if 



exp 



cos 



V - p) 1 

(p 2 - 



i] 



for v > p 

for p > v 

for v > p 

^ + 1] for p> V . 



- (42) 



{p — v ) 



1x1 «1- 



(43) 



Introducing these approximations for the Bessel functions in both 
equations (35) and eliminating \p between them, we obtain the char- 
acteristic equation for the bent slab 



Sp 



h t(p?» - " 2 ) i - «• - " 2 ) § ] 



= pir — tan" 



Wi 



P13 — v 

2 2 

L" — P3 J 



1 + i exp 



2> 2 - 
3 



**$ 



.v, 2 / 2 2 \i 

-tair'^lfe— r 
^i \" ~~ pi 



(44) 
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in which p is an arbitrary integer bigger than zero which determines 
the order of the mode in the x direction, and the tan -1 functions are 
to be taken in the first quadrant. 

Let us rewrite this equation substituting P3 , P5 , p i3 , and pis by the 
values given in equation (36) ; furthermore, let 



A* = 



and 



v = k t R 
K = (fc? - K - &. 



(45) 
(46) 

(47) 



Because of these two last definitions, k„,k a , and k y are the axial and 
the transverse propagation constants at x = 0. The characteristic equa- 
tion (44) then becomes 



Rk 
3k 



»[■-( 



- 1 - 



klR) J 
1 + t 



= pir — tan 



_. n* 



j 2 t 3 R [\ (k x A 3 \ a ~\ 



tan —3 

n, 



1 - 



B. 



(*; - *s - « 



«-(i-|)«-«) 



(48) 



To solve this equation for k a we expand the left side and the second 
tan -1 in powers of 1/R and the first tan -1 in powers of the exponential. 
Assuming R is large and keeping the first term of each perturbation 
calculation, the solution of equation (48) is 



, | , , 2c .k l0 a, 



where 



c = 



7TO 



1 1 + 2F S 



2fc l0 a U,/ (R 1 + F 3 + Fs 



(49) 



(50) 
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and 



fczll 



1 - 



k,«A. 



e v p 



it-w-y( 



ak x J J 



1 + ^3 + F* 



in which 



^ = 



A : 



TTCl 



1 - 



k x0 A< 



1 - 



2tt 3 K 



1 g 

A 2 

k 3 R 



k x0 Ai 



<R == firTi = 2(ni - nj) 1 . , 

ft/, o-^3 rc »0 

"'JO = ("'l "'II "'io) 1 

and fc.ro is the solution of the equation 



(51) 
, (52) 

(53) 

(54) 



n 

k z0 a = pw — tan -j r 

7t 



1 



, -in, 
— tan -s t= 






(55) 

This is the physical interpretation of equation (49) : the trans- 
verse propagation constant k x measured at x = is made of three 
terms. The first term, kjo, is the transverse propagation constant of 
the guide without curvature ; the second and third terms are perturba- 
tions related to the change of field profile and radiation introduced by 
the curvature. It is easy to find that c" is the mode conversion loss that 
would exist at a junction between a straight guide and a curved one, 
and a c is the attenuation coefficient of the curved guide. 

The field components in media 3 and 5 decay almost exponentially 
away from the guide. The length $ s , over which the intensity in medium 
3 decays by 1/e, is derived as in equation (41) to be 



£3 = 
and only approximately 



k Z 3 



\k\ k 3 I k x 



(56) 



fc - 



In 



(k;-k 2 5 - I AS|)»' 



(57) 
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All these equations have been derived under the assumption that 
inequality (43) is satisfied; this means that the field configuration of 
the curved guide is very close to that of the straight guide. In other 
words, c <SC 1. For a given R, if one chooses the width a of the guide 
large enough, these inequalities are not satisfied, the previous results 
are no longer applicable, and a new solution is needed. We proceed 
to find it. 

Let us assume as a limiting case that in Fig. 2 

a = R. (58) 

The characteristic equation derived from the first equation of (35) , 
making \p = 0, is 

(pi3 - v 2 ) 1 . i\ +«,-» ?1 sY pi3 ~ " 2 Y 
KH \ 2 — = (P — t)t - tan -a I— ~2 ) 



\ l+iew [-\^q. (59 ) 



Following similar steps to those taken to solve equation (44), we 
substitute p 16 , p a , and v by the values given in equations (36) and 
(46) ; we obtain 



w? = (v " * )7r 



— tan 



-1 ^3 



1 +t 



1 exp nggrlTJ. 



[fer - ■]' 



(60) 



The primes distinguish the symbols from those used previously. 
To solve this equation we notice that for small losses it must be that 



fcx-Ai 



«1. 



(61) 



Therefore, the tan -1 can be replaced by its argument and the approxi- 
mate solution of equation (60) is 



K = K 



(icy} ' 



(62) 



where 



W3 K 1 q 

£J kR{n\ - nf)* 
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•exp 



CR' 



1 - 



6r(p - fi 
(ft' 



2n' 



1 _ - 1T2 
1 „ ■} 



6 



3 n\ W(p - i) 2 rti 



k', = [*■ - « - (#„)*]», 

* ~.4 3 |_ (R' .l 1_ 3flfL? 



(p " \Y<K- 



and 



(R = 777T2-T3 = 2(n, - n 3 )' jyr- 



QUrAi 



(KoY 



The field components outside the guide decay to 1/e in a length 
1 1 



S- 



K, 



I \2-il ' 



VK -K- (A-; o ) 2 ] 



(63) 

(64) 

(65) 



(66) 



(67) 



a. 2 E" VQ Modes: Polarization Along y 

The field components and propagation constants can be derived from 
those in Section A.l by changing E into H, n into — e, and vice versa. 
Except for their polarizations, the El Q and E" vq modes are very similar. 

The formulas equivalent to equations (37) and (41) are 



K'b = qir — tan ' -§ j= 



AJP> ~ ' 



-, - tan -g F 
1 n 



_\^4 4 /c l 



- 1 



(*? - « - (fc-) 2 ) 1 



(68) 
(69) 



The double prime distinguish these symbols from those used before. 

The equivalent formula to any of those between equation (45) and 
(67) can be derived from that formula by substituting the ratio of 
refractive indexes by unity, but leaving the differences between squares 
of indexes unchanged. For example, the formula equivalent to equation 
(52) for E" modes is 



Fi' = 



; 


Ic" A 
" j- u ■< i 3 


2~ 


5 


ira 1 — 

L 


5 

7T 







1 - 



w 3 

1 - -I 



fcxoAi 



(70) 
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